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In this paper we examine whether indirect detection constraints on dark matter associated with a non-thermal 
history may be significantly improved when accounting for the presence of galactic substructure in the form of 
dark matter spikes. We find that significant constraints may be derived from the non-observation of an excess of 
diffuse gamma-rays and from the properties of bright gamma- ray point sources observed by the Fermi Gamma- 
Ray Space Telescope, but these constraints depend sensitively on the details of the formation of the first stars 
and their subsequent black hole remnants. However, we also find that, especially if WIMPs annihilate primarily 
to quarks or gauge bosons, it is possible to extract meaningful and conservative bounds on the annihilation cross 
section. 
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I. INTRODUCTION 

Despite the successes of precision cosmology, ex- 
isting observations seem to tell us little about the his- 
tory of the universe prior to Big Bang Nucleosynthe- 
sis (BBN). This is unfortunate given the expectations 
from particle theory for a rich amount of phenomenol- 
ogy at these scales; including symmetry breaking tran- 
sitions, the generation of mass in the Standard Model, 
the origin of the baryon asymmetry, and the existence 
of Cold Dark Matter (CDM). The Large Hadron Col- 
lider (LHC) is cuiTently probing the microphysics re- 
sponsible for many of these processes, but without an 
exact description of the cosmological history prior to 
BBN our understanding will remain incomplete. 

In the case of CDM, the standard assumption of a 
thermal history prior to BBN provides a well-motived 
and convincing scenario for connecting the cosmolog- 
ical and microscopic origin of CDM [1]. In this ap- 
proach, one assumes that very early in its history the 
universe achieves thermal equilibrium and remains in 
that state until the time of BBN. In such an approach, 
the amount of CDM today depends parametrically on 
the properties of the CDM particles (mass and cross 
section) and the temperature at which the particles 
ceased to annihilate - so-called 'freeze-out'. It is re- 



assuring that when comparing this estimate with pre- 
cision cosmological measurements for the amount of 
CDM today we get a prediction for the mass and an- 
nihilation cross section near the scale of electroweak 
symmetry breaking. However, in the simplest models 
this "WIMP Miracle" is spoiled by a tension with elec- 
troweak precision constraints. Given that search strate- 
gies at LHC and other CDM detection experiments de- 
pend on assumptions about the self-annihilation cross 
section of these particles' it is crucial to establish the 
robustness of the thermal scenario and the associated 
cosmological constraint, i.e. {(Jv)th ~ lO^'^® cm'^s"'^, 
as well as identifying any other viable alternatives for 
their production [2|1 ■ 

One possible alternative scenario is that of a non- 
thermal history. This scenario occurs if massive par- 
ticle decays or phase transitions lead to a significant 
entropy and particle production prior to BBN. If such 
transitions occur after the thermal freeze-out of CDM, 
predictions for the microscopic properties of the total 
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' There is both a direct connection, as is the case for the indirect 
detection of annihilation products, as well as a more implicit con- 
nection that appears in model-specific scenarios. As an example of 
the latter, in the Minimal Supersymmetric extension of the Stan- 
dard Model (MSSM), requiring a thermal dark matter candidate 
typically leads to a WIMP that is a bino-like neutralino. However, 
if one drops the thermal constraint more regions of the MSSM pa- 
rameter space become viable, which can lead to different possible 
benchmark signatures at LHC |[^. 
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amount of CDM may differ significantly from the usual 
thermal scenario IH]. Such scenarios have deservingly 
received much skepticism over the years, particularly 
because the non-thermal production of CDM must oc- 
cur in a very narrow window - after CDM thermal 
freeze-out but prior to the onset of BBN - naively intro- 
ducing a new and unmotivated scale of physics into the 
problem. However, in the particular case of Anomaly 
Mediated Supersymmetry (SUSY) breaking one finds 
that this scale is set by the scale of SUSY breaking 
and non-thermal CDM is a natural prediction for this 
class of models |4]. Building on this intuition, more 
recently it has been suggested that this may be a gen- 
eral expectation of a larger class of gravity-mediated 
SUSY models when one accounts for theoretical self- 
consistency in the ultraviolet ||3ti5t]. In addition to this 
theoretical motivation, non-thermal models make def- 
inite and testable predictions which are currently be- 
ing scrutinized at colliders, as well as by ground- and 
space-based CDM searches. One such prediction is the 
enhancement of the self-annihilation cross section by 
as much as three orders of magnitude compared with 
that of the standard thermal scenario, while still yield- 
ing the correct amount of CDM cosmologically ^ . 

In this paper we focus on whether indirect detection 
constraints on non-thermal CDM can be significantly 
improved when accounting for galactic substructure re- 
sulting from dark matter spikes. Specifically, we fol- 
low in their analysis of the gamma-ray constraints 
on dark matter annihilation in the dark matter spikes 
in our Galactic halo. Dark matter spikes arise due to 
the contraction of a dark matter minihalo when a bary- 
onic object (e.g. a star) forms at its center, as was the 
case with the first generation of stars to form in our uni- 
verse. Indeed, as a result of the increased dark matter 
density in the spike, the very first stars are thought to 
have undergone a phase during which they were sup- 
ported by dark matter annihilations, dubbed the Dark 
Star phase [7]. The affect of a boosted annihilation 
cross section on the evolution of Dark Stars was exam- 
ined in ll^, where they found that the Dark Star phase 
is shortened by an enhanced dark matter annihilation 
cross section, though the existence of the phase is ro- 
bust. 



PAMELA^ excess is in-fact a signature of dark matter 
annihilations, the data suggests that dark matter anni- 
hilates preferentially to leptonic final states [J3i1 . It is 
possible to construct such models lfT2ll . however it is 
noteworthy that predicted fluxes of charged particles 
can suffer from large uncertainties associated with as- 
tro|)hysical backgrounds. Indeed, it was demonstrated 
in [91] (see also 1 14,15"]) that in the case of anti-protons, 
astrophysical backgrounds can be significantly lower 
than previously expected while still being consistent 
with the Boron to Carbon ratio. 

Here we examine the potential of the gamma-ray 
data from the Fermi Gamma-Ray Space Telescope 
(FGST) to constrain models of non-thermal dark mat- 
ter in local spikes. Uncertainties in the astrophysical 
backgrounds play an inconsequential role in the fol- 
lowing analysis, however, we find that constraints on 
dark matter annihilation can be ambiguous in the ab- 
sence of a reliable star formation history. We empha- 
size that an important and difficult challenge for this 
program is the establishment of constraints on the typ- 
ical mass and formation era of the first generation of 
stars. For example, upcoming observations with the 
James Webb Space Telescope (JWST) may provide 
some hints about the formation of the first stars lfT6l . 
and it's even possible that JWST will observe a Dark 
Star 111?]. As we demonstrate, if a star formation his- 
tory is established, the constraints on dark matter anni- 
hilation in local spikes may be very significant. 

In the next section we briefly review the mechanism 
by which dark matter spikes form in the early universe 
and the method used to extract the local distribution 
of surviving spikes in our Galactic halo. In section Hill 
we demonstrate how the presence of this substructure 
can lead to a general enhancement of the gamma-ray 
constraints on non-thermal dark matter model building. 
We briefly conclude in the section that follows. 



11. GALACTIC SUBSTRUCTURE FROM CDM 
SPIKES 

One might expect that the distribution of CDM 
within our galaxy can be strongly influenced by the for- 
mation and evolution of objects such as black holes. In 



Current constraints from indirect detection already 
put strong bounds on the allowable cross sections 
for non-thermal models JH [HI]. Additionally, if the 
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particular, as the gravitational potential becomes dom- 
inated by a compact baryonic object, the CDM distri- 
bution near this object will be affected. Gondolo and 
Silk have examined this possibility for the supermas- 
sive back hole at the center of our galaxy (around 10^ 
M0) around which one might expect a large enhance- 
ment in the CDM density lUSll . However, further inves- 
tigations revealed that such extreme inhomogeneities 
are most likely negligible today due to a number of ef- 
fects, including major merger events, off-center forma- 
tion of the seed black hole, gravitational scattering off 
stars, and CDM annihilations 1 19-22]. Zhao and Silk 
then proposed |23] that these wash-out effects may not 
be present for small over-densities, or spikes, resulting 
from Intermediate Mass Black Holes (IMBHs), which 
are the expected remnants of the earliest stars to form, 
known as Pop-Ill stars^. 

Bertone, Zentner, and Silk (BZS) examined this pos- 
sibility in more detail in f24] (for a review see ||25| ) 
using an analytic model of halo evolution and perform- 
ing 200 statistical realizations for the growth of a Milky 
Way-sized halo. The population of IMBHs was gener- 
ated by identifying Scr over-densities in the smoothed 
primordial density field at a redshift of z = 18 and re- 
placing each of those peaks with a 100 Mq black hole. 
Tracking the growth and mergers of the structures un- 
til today, they find an expected number of IMBHs in 
our galaxy to be Nbh = 1027 ± 84. The uncertainty 
in this number reflects unknowns in the model parame- 
ters, such as the redshift at which small-scale fragmen- 
tation of baryonic disks becomes important and black 
hole seeds cease to form. BZS accounted for this un- 
certainty by varying the redshift at which the seeds are 
initially evolved. 

A different approach was taken by one of us (PS) in 
collaboration with J. Diemand, K. Freese, and D. Spol- 
yar in (see also [26.] ). hereafter referred to as SDFS. 
Their analysis uses the Via Lactea-II cosmological N- 
body simulation \2f] to estimate the number and mass 
distribution of CDM minihalos as a function of red- 
shift. Minihalos suitable for the formation of Pop-Ill 
stars are identified at high redshift, and the distribution 
evolved until today. Assuming each minihalo hosted a 
Pop-Ill star, these minihalos exist today as spikes, each 
surrounding an IMBH Pop-Ill remnant. The distribu- 
tion of these IMBHs and surrounding spikes depends 



^ In this paper we will carelessly refer to Population III.l as Pop-Ill. 



on the duration of Pop-Ill star formation, though the 
exact redshift at which Pop-Ill star formation ceases 
remains uncertain. Additionally, the density profile of 
each individual spike, and therefore the expected dark 
matter annihilation rate, depends on the typical size of 
the remnant black holes. Here we build on the analysis 
of SDFS, whose methodology we now review, refer- 
ring to the original papers for more detail 

Via Lactea-II ]27] is the first cosmological N-body 
simulation of a MiUcy-Way sized dark matter halo capa- 
ble of resolving the ^10^ Mq minihalos in which the 
first stars formed. Star formation depends on the ability 
of the baryonic clouds to efficiently cool as they col- 
lapse. This cooling proceeds primarily through excita- 
tions of molecular hydrogen, the abundance of which 
depends on the temperature and therefore redshift of 
formation. Using this fact, Trenti and Stiavelli found 
a minimal mass for minihalos in which Pop-Ill stars 
could have formed of ll28|] 

/, X -2.074 

M^i: = 1.54 X 10^ Mq (^^j . (1) 

The maximum mass of halos that formed Pop-Ill stars 
is less important since the hierarchical nature of struc- 
ture formation favors small mass minihalos, but for 
completeness SDFS took a maximum mass of 10'' Mq. 
Given the uncertainty in the redshift at which Pop- 
III star formation gave way to the formation of less- 
massive subsequent generations of stars (which are 
not expected to result in the spikes we examine here), 
SDFS considered three possible termination redshifts 
Zf = 11, 15, and 23. For brevity, here we consider 
only Zf — 15. 

Assuming each Pop-Ill star ended its life by collaps- 
ing to a black hole, and given a Pop-Ill termination 
redshift and the viable minihalo mass range above, the 
current number density of black holes surrounded by 
spikes, Nbfi, is related to the total possible number of 
viable minihalos, Nhaios, by 

^6/1 = /o (1 ^ fmerged) Nfialos, (2) 

where /o is the fraction of halos that are expected 
to host Pop-Ill stars, and fmerged is the fraction of 
CDM spikes that are destroyed by black hole mergers. 
SDFS argued that mergers are most important for the 
highest mass black holes and for /o ~ 1, in which 
case they would reduce the number of spikes by at 
most a factor of two. For lighter black holes and/or 
smaller /o, it was argued that this effect is negligible 
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and N, 



bh 



foNhal. 



Here we will fix the fraction of 



black holes to form and survive, fg — /o(l — f merged) 
and consider two possible values; fg — 0.1 and the 
maximal case = 1. 

If the growth rate of a baryonic object at the center 
of a minihalo is slow with respect to the time it takes 
CDM particles to cross the central region, the contrac- 
tion of particle orbits and the formation of CDM spikes 
may be modeled by adiabatic contraction. SDFS used 
the Blumenthal et al. prescription for adiabatic con- 
traction fl^, which predicts a roughly power-law den- 
sity profile for the spikes and is independent of the 
CDM particle mass. However, given the enhanced den- 
sity of CDM, some particle self-annihilations will take 
place. This depends on the lifetime of the central mass 
core (tbh) and leads to an upper limit on the CDM den- 
sity 



{(Jv)tbh ' 



(3) 



where and (av) are the CDM mass and averaged 
self-annihilation cross section times velocity, respec- 
tively. 

In summary, SDFS find that the formation of Pop- 
III stars leads to a significant number of CDM spikes 
in our own galaxy today, as first anticipated by Zhao 
and Silk. In the next section, we consider the feasi- 
bility of using the existence of these spikes to sharpen 
constraints on the properties of CDM through (non- 
)observation of their annihilation products by FGST. 



III. NON-THERMAL CDM CONSTRAINTS FROM 
SPIKES 

For a Majorana CDM particle with mass and 
average annihilation cross section times velocity {(Jv), 
the rate of self annihilations inside a spike is 



{av)_ 
2rn^ 



dr 47rr2 p%,k^{r), 



(4) 



where Tmin and r.max are the inner and outer radii of 
the CDM spike in which annihilations occur with the 
former being of order the Schwarzschild radius of the 
black hole, and Pspike {r) is the CDM density profile of 
the spike. 

We consider several WIMP candidates defined by 
their masses and annihilation channels. Calculations 
are performed for WIMP masses of 100, 200, 500, 



1000, and 2000 GeV and Standard Model final states 
bb, W~^W~, T^T^, and p^ p^ ■ The resulting spec- 
trum of photons dNf/dE from annihilation to final 
state / is computed with PYTHIA |30]. For xx ~^ 
the photon spectrum comes from final state ra- 
diation and is given by [3 1.1 



dN„ 



dx 



- 2x 



xn/a 



In 



s{l-x) 



(5) 

where x = E^/m^, the center-of-mass energy squared 
is s = 4to^, and a « 1/137 is the fine structure con- 
stant. We note that WIMP candidates typically annihi- 
late to a variety of final states with the rate of annihila- 
tions in a CDM spike expressed as 



Tf = BfV 



(6) 



where i?/ is the branching ratio to the final state /. The 
intrinsic photon luminosity from CDM annihilations in 
any CDM spike is then 



C= [ dEj2 

f 



dNf 
dE 



(7) 



Given the luminosity resulting from dark matter annihi- 
lations in spikes, we can now proceed to establish con- 
straints on the CDM self annihilation cross section us- 
ing both point source and diffuse flux data from FGST. 



A. Point Source Constraints 

We first consider establishing constraints on the 
WIMP self annihilation cross section by requiring that 
annihilations in the nearest spike do not lead to a 
point source flux that exceeds that from the bright- 
est recorded FGST point source. Point source con- 
straints rely heavily on the estimate of the distance to 
the nearest spike, determined by integrating the prob- 
ability density of finding a spike in the neighborhood 
of our Solar System. Despite the fact that the bright- 
est FGST point source is associated with the Vela pul- 
sar l32ll . in this analysis we simply require that the 
gamma-ray flux from the brightest spike not exceed the 
gamma-ray flux from Vela, resulting in the somewhat 
bizarre requirement that the brightest spike must be lo- 
cated along our line-of-sight to Vela. As this possibility 
is not excluded, we reserve further discussion of this is- 
sue until the end of the section. 
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As discussed in section [III given the uncertainties 
in spike formation we will consider both small mass 
{msH = 100 M0) and large mass (mBH = 10"* Mq) 
black holes and we will consider two values for the 
fraction of black holes to form and survive; fg ~ 0.1 
and /s = 1. Our results for the point source analy- 
sis are given in Figures [T] and |2] In Figure [T] we fix 
fs = 0.1 and present the upper limit on the average 
WIMP annihilation cross section times velocity as a 
function of the mass of the dark matter particle for a 
typical black hole mass of 100 M© (top panel) and 
10"* Mq (bottom panel) for four choices of final state 
particles; bb (solid black curves), W^W^ (dashed 
black curves), jj.'^ (solid grey curves), and r+r^ 
(dashed grey curves). In each panel we present the cos- 
mologically determined thermal WIMP cross section 
{<^v)th = 3 X 10^^^ cm^s^^ for comparison. Simi- 
larly, in Figure 12] we show the upper limit on {av) as a 
function of WIMP mass for fs = l. 

It is clear that the constraints are quite sensitive to 
both the dark matter annihilation mode and the typical 
black hole mass. Constraints are also sensitive to the 
value of fs, an as-yet unknown parameter The least 
constrained case considered here is shown in the top 
panel of FigurefT] with rribh — 100 Mq and fs — 0.1. 
Comparison of these limits with the left panel of Fig- 
ure 5 in i33i1 reveals that the constraints are compa- 
rable^. However, for the heaviest black holes consid- 
ered here and fs — 1, as shown in the lower panel of 
Figure |2] we see that non-thermal models with light 
WIMPs are essentially completely ruled out, with the 
possible exception of the case where XX ^ M^A* ■ 
We expect that the true effect of the spikes will lie 
somewhere between these two extreme cases. How- 
ever, we see from the upper panel of Figure [T]that even 
in this case meaningful constraints may be achieved for 
WIMPs with final states of bb and W+W'. The lat- 
ter places an important constraint on the non-thermal 
wino-like LSP scenario discussed in ^ [l^, which 
used the non-thermal enhancement of the cross section 
to address the PAMELA data. 

Finally, it is obvious that the constraints presented 



' Differences in the slope of the constraint in the (m-^ , (crv) ) plane 
are due in part to the fact that the spike profiles (specifically Pmax ) 
are affected by the WIMP mass. As WIMP mass increases, not 
only does the luminosity of a particular spike decrease due to the 
decreased number density of WIMPs, but the luminosity also de- 
creases because the spike has essentially lost more of its core. 
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FIG. 1: Upper limit on the annihilation cross section as a 
function of the mass of the dark matter particle for a typical 
black hole mass of 10^ Mq (top panel) and 10^ Mq (bot- 
tom panel) for four choices of final state particles; bb (solid 
black curves), W^W~ (dashed black curves), fi'^ p~ (solid 
grey curves), and t~^t~ (dashed grey curves). Here we as- 
sume fs ~ 0.1. The horizontal line in each panel indicates 
{av)th = 3 X 10~^* cm^s'^. 



here would improve had we chosen to use the flux limit 
from the brightest unassociated FGST point source, 
rather than the flux limit from Vela. For a spike 
located at some distance D from our solar system, 
$ oc {av)/D'^, so it is possible to translate the con- 
straints on {(jv) from Vela to constraints from the 
brightest unassociated source, which has an integrated 
luminosity ^ 1/22 that of Vela. Therefore the limits 
from the brightest unassociated source, for the choices 
of fs in Figures [T] and |2] would simply be shifted to 
lower {(jv) by a factor of ~ 1/22, representing a 
notable improvement in the ability to constrain non- 
thermal cross sections. Alternatively, one could imag- 
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FIG. 2: Upper limit on the annihilation cross section as a 
function of the mass of the dark matter particle for a typical 
black hole mass of 10^ Mq (top panel) and 10* Mq (bottom 
panel) for four choices of final state particles; bb (solid black 
curves), W^W' (dashed black curves), fi^ fJ.~ (solid grey 
cun'es), and t^t~ (dashed grey curves). Here we assume 
fs ~ 1. The horizontal line in each panel indicates (av)th ~ 



ine that fs is in fact less than the minimal value of 0.1 
that we have chosen to examine here. If we derive our 
limits according to the flux from the brightest unasso- 
ciated FGST point source, the constraints in Figures [1] 
and|2]would apply to fs ~ 0.001 rather than 0.1 (upper 
panels), and fs w 0.01 rather than 1 (lower panels). 
If it is true that there are no bright spikes located along 
our line of sight to any of the brightest associated FGST 
point sources, then even if fs is quite small there are 
very significant limits on (av) from non-observation 
of bright nearby CDM spikes. 

With these examples of how point source flux from 
WIMP annihilations in spikes can be used to place con- 



straints on non-thermal dark matter, we now turn to the 
case of the diffuse flux. 



B. Diffuse Flux Constraints 
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FIG. 3: Cross sections that may be excluded by the FGST 
measurement of the diffuse gamma-ray flux for fs — I and 
annihilations to fj,'^ n~ . If the typical size for a black hole 
is 100 Mq, the shaded region between the solid contours is 
excluded. If the typical size for a black hole is 10* Mq, the 
region between the dashed contours is excluded. The hor- 
izontal line indicates the thermal dark matter cross section 
{av)th = 3 X 10"^*^ cm^s-^. 

The CDM annihilation cross section may also be 
constrained by requiring that the diffuse flux from dark 
matter annihilations in the spikes not exceed the FGST- 
measured diffuse gamma-ray flux by more than 3cr in 
any of the nine energy bins of Ref. |34]. This pro- 
vides a quite robust constraint given our adoption of the 
conservative assumptions for the number of spikes that 
contribute to the diffuse flux as outlined in Ref. |26]. In 
Figure |3] we present the values of (cru) that may be ex- 
cluded by the diffuse gamma-ray flux as measured by 
FGST for the channel XX n+fj,- for fs ^ I. The 
shaded region represents the cross sections that are ex- 
cluded if the typical black hole mass is 100 M0, while 
the region between the dashed contours is excluded if 
the typical black hole mass is 10** M©. As expected, 
cross sections that are below the accessible range result 
in too low a photon flux to provide a meaningful con- 
straint. However, cross sections above the accessible 
range result in extremely bright spikes, such that many 
or most of the spikes in our Galactic halo would be 
visible as point sources, and therefore very few would 
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contribute to the diffuse flux. The horizontal line repre- 
sents the standard cross section for thermal dark matter, 
{o'v)th = 3 X lO^'^^ cm'^s"^. We see that even for the 
case of small black hole mass that light non-thermal 
WIMPs with the XX f^~^IJ'~ channel may lead to 
significant constraints. We remind the reader, how- 
ever, that the largest diffuse flux is expected in models 
where the luminosity of an individual spike is very low 
(thus many/most spikes contribute to the diffuse flux). 
Therefore, of the cases considered here, XX M^M" 
with /s = 1 results in the strongest diffuse constraints. 
For fs < 1, the diffuse constraint weakens. As was 
concluded in |6], there are few cases where the dif- 
fuse constraint is stronger than that from point source 
brightness. 

IV. CONCLUSIONS 

We have examined whether significant constraints 
on non-thermal dark matter can be derived by account- 
ing for the presence of dark matter spikes in our Galac- 
tic halo. We find that, despite the uncertainties in 
the formation of dark matter spikes and the associated 
black holes, meaningful constraints may be expected 
for fs ^ 0.1, and even for fs as small as 10^"^ if 
one is willing to accept that there are no spikes hid- 
ing along our line of sight to the brightest associated 
gamma-ray point sources. Existing constraints on non- 
thermal dark matter annihilation cross sections may be 
improved by the non-observation of a gamma-ray flux 
from spikes in our Galactic halo, especially if the typi- 



cal mass of a black hole at the center of a spike is rather 
large 10^ M©). We have also demonstrated that 
constraints can be established based on the contribu- 
tion of faint CDM spikes to the diffuse gamma-ray flux 
for the example case of xx ~^ A'^M ■ For this particu- 
lar final state, and if each minihalo capable of forming 
a Pop-Ill star did form one, we find that non-thermal 
WIMPs are restricted to be quite massive, even for the 
lighter 100 M0 central black holes. 

Despite the many uncertainties in the star formation 
history, these results are promising and merit further 
investigation into the importance of dark matter spikes 
in indirect detection of non-thermal dark matter How- 
ever, an important and difficult challenge for this pro- 
gram is to better establish the typical mass and forma- 
tion era of the first generation of stars. As we have 
shown, if a star formation history is better established, 
accounting for dark matter annihilation in local spikes 
may significantly improve existing constraints on non- 
thermal dark matter 
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